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ABSTRACT 

The Angstrom Project is using a global network of 2m-class telescopes to conduct a high cadence pixel 
microlensing survey of the bulge of the Andromeda Galaxy (M31), with the primary aim of constraining its 
underlying bulge mass distribution and stellar mass function. Here we investigate the feasibility of using 
such a survey to detect planets in M31. We estimate the efficiency of detecting signals for events induced 
by planetary systems as a function of planet/star mass ratio and separation, source type and background M3 1 
surface brightness. We find that for planets of a Jupiter-mass or above that are within the lensing zone (~ 
1-3 AU) detection is possible above 3 a, with detection efficiencies ~ 3% for events associated with giant 
stars, which are the typical source stars of pixel-lensing surveys. A dramatic improvement in the efficiency of 
~ 40 - 60% is expected if follow-up observations on an 8m telescope are made possible by a real-time alert 
system. 

Subject headings: gravitational lensing - planets and satellite: general - galaxies: individual (M31) 
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1. INTRODUCTION 

Various techniques are being used to search for ex- 
trasolar planets, including the radial velo city te chnique 
JMavor & Ouelozl 119951 iMarcv & Butted 119961) tran- 
sit method (Struve direct imaging llAngell 1 1994: 

Statu & Sandlerl 119951) . pulsar timing JWolszczan & Fraill 
19921). an d mic rolensing ilMao & P aczYflskj I199J1 
Gould & Loeblll992l) . See the reviews of lPerrvmad J20001 



2005). The microlensing signal of a planetary companion 
to microlens stars is a short-duration perturbation to the 
smooth standard light curve of the primary-induced lensing 
event occurring on a background source star. Once the signal 
is detected and analyzed, it is possible to determine the 
planet/star mass ratio, q, and the projected planet-star sepa- 
ration, s (normalized by the angular Einstein ring radius 6 E ). 
Recently, two r obust microl ensing det ections of exo planets 
were reported bv lBond etal] J2004I) and lUdalski et alJ ( l2005T) . 

The microlensing technique has various advantages over 
other methods. First, microlensing is sensitive to lower- 
mass planets than most other methods (except pulsar timing) 
and it is possible, in principle, to det ect Earth-mass plan- 
ets fr om ground-based observations ( Gou ld. Gaudi. & Hani 
2004). Second, the microlensing technique is most sensitive 
to planets located in the so-called lensing zone corresponding 
to the range of 0.6 - 1.6 Einstein ring radii. The typical value 
of the Einstein radius, rg, for Galactic lensing events is a cou- 
ple of AU, and thus the lensing zone roughly overlaps with the 
habitable zone. Third, the microlensing technique is the only 
proposed method that can detect and characterize free-floating 
planets jBennett & Rhiell2002t lHan et alJl2005l) . Fourth, the 
biases in the search technique are less severe and can be quan- 
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tified easily compared to other methods JGaudi et alJl2002|) . 
Therefore, the microlensing technique will be able to provide 
the best statistics of the Galactic population of planets. 

In addition to the advantages mentioned above, the mi- 
crolensing technique is distinguished from other techniques 
in the sense that the planets to which it is sensitive are much 
more distant than those found with other techniques. With 
the advent o f photometry techniques like the diffe rence imag- 
ing (jAlard & Luptonl 19981) and pixel method ( Melchi or et alJ 
1999), microlensing searches are not restricted to the field 
within the Galaxy and can be extended to unresolved star 
fields of nearby galaxies such as M3 1 . Therefore, microlens- 
ing is the only feasible technique that can detect planets lo- 
cated in other galaxies. 

Microlensing searches toward M3 1 have been and are being 
carried out by various co l labora ti ons including t he POINT- 
AGAP E (lAuriere et alJ 1200 It iPaulin-Hendriksson et alJ 
120021 I2003t lAneTatl 12004 Belokurov et alJ 120051) . 
AGA PE JAnsarietalJ [19971 1 19991). VATT -Colombia 
llCrotts & Torrmnevl 119961 lU glesich et aU 120041). MEGA 
ide Jong et alJ 120041) . and WeCAPP jRiffeser et all 120011 



2003) collaborations, as well as MDM (Calchi Novati, et al 
2003), McGraw-Hill llCalchi Novati. et alJ 120021) . and 
Nainital (Uoshiet all 120031120051) ' surveys. The monitoring 
frequencies of these experiments are typically ~ 3 observa- 
tions per week, too low to detect planetary signals. However, 
with the expansion of the global telescope network, the mon- 
itoring frequency of such surveys is rapidly increasing. For 
example, a new M31 pixel-lensing survey, the Andromeda 
Galaxy Stellar Robotic Microlensing (Angstrom) project is 
expected to achieve a monitoring frequency of ~ 5 observa- 
tions per 24-hour period by using a network of telescopes, 
including the robotic 2m Liverpool Telescope at La Palma, 
Faulkes Telescope North in Hawaii, 1.8m telescope at the 
Bohyunsan Observatory in Korea, and the 2.4m H iltner 
Teles cope at the MDM Observatory in Arizona (iKerins et al j 
2005j). 

The possibility of detecting planetary microlensing events 
caused by a lens loc ated in M31 was discussed by 
iBaltz & Gondol<3 (1200 ll) . However, the main focus of that 
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paper was evaluating the detectability of events caused by bi- 
nary lenses in general and the comment about the planetary 
lensing was brief treating the planetary system as one case of 
the binary lenses. In addition, their detection rate estimate of 
the M31 planetary lensing events was based only on events 
that exhibit caustic crossings, while a significant fraction of 
events with detectable planetary signals might be non-caustic- 
crossing events. Moreover, the aim of their work is was the 
rough evaluation of feasibility and thus not based on specific 
observational setup and instruments. Similarly, the work of 
ICovone et al.l (120001) was also based on an arbitrary observa- 
tional setup. 

In this paper, we explore the feasibility of detecting planets 
in M3 1 from a high-frequency pixel-lensing survey using a 
global network of 2m-class telescopes. The paper is organized 
as follows. In § 2, we briefly describe the basics of planetary 
lensing. In § 3, we estimate the efficiency of detecting plane- 
tary signals for events induced by planetary systems with var- 
ious planet-star separations and mass ratios, associated with 
source stars of different types, and occurring toward fields 
with a range of surface brightness //. From the dependence 
of the detection efficiency on these parameters, we investigate 
possible types of detectable planets, the optimal source stars, 
fields, and observation strategy for M3 1 planet detections. In 
§ 4, we discuss about methods to further improve the planet 
detection efficiency. We summarize the results and conclude 
in §5. 

2. BASICS OF PLANETARY LENSING 

Planetary lensing is described by the formalism of a binary 
lens with a very low-mass companion. Because of the very 
small mass ratio, the planetary lensing behavior is well de- 
scribed by that of a single lens of the primary star for most 
of the event duration. However, a short-duration perturbation 
can occur when the source star passes the region near a caus- 
tic, which represent the set of source positions at which the 
magnification of a point source becomes infinite. The caus- 
tics of binary lensing form a single or multiple closed figures 
where each figure is composed of concave curves (fold caus- 
tics) that meet at cusps. 

For the planetary case, there exist two sets of disconnected 
caustics. One 'central caustic' is located close to the host star. 
The other 'planetary caustic' is located away from the host 
star and its number is one or two depending on whether the 
planet lies outside (s > 1) or inside (s < 1) the Einstein ring. 
The size of the caustic, which is directly proportional to the 
planet detection efficiency, is maximized when the planet is 
located in the 'lensing zone', which represents the rang e of th e 
star-planet separation of 0.6 < s < 1.6 (Gould & Loeb 1992). 

The planetary perturbation induced by the central caustic 
is of special interest for M31 pixel-lensing events. While the 
perturbation induced by the planetary caustic can occur at any 
part of the light curve of any event, even those of low magnifi- 
cation, the perturbation induced by the central caustic always 
occurs near the peak of the light curve of a high-magnification 
event. Then, the chance for the M31 pixel-lensing events 
to be perturbed by the central caustic can be high because 
these events tend to have high magnifications. In addition, 
the chance of detecting planetary signals for these events be- 
comes even higher due to the improved photometric precision 
thanks to the enhanced brightness of the lensed source star 
during the time of perturbation. 

3. DETECTION EFFICIENCY 



To estimate the efficiency of detecting planetary signals of 
M3 1 events, we compute the 'detectability' defined as the ra- 
tio of the planetary signal, e, to the photometric precision, cr p h, 



i.e., 



V=±±. 



(1) 



The planetary signal is the deviation of the lensing light curve 
from that of the single lensing event of the primary lens star, 
and thus it is defined as 

A-A 



e = 



(2) 



where A is the magnification of the planetary lensing and Aq is 
the single lensing magnification caused by the host star alone. 
For an M3 1 pixel-lensing event, the lensing signal is the flux 
variation measured on the subtracted image, while the noise 
is dominated by the background flux. Then, the photometric 
precision can be approximated as 



(3) 



F s (A-iy 

where F$ and Fb are the baseline flux of the lensed source 
star and the blended background flux, respectively. Under this 
definition of the detectability, V = 1 implies that the planetary 
signal is equivalent to the photometric precision. 

We estimate the detection efficiency for a representative 
event that is most probable under the assumption that the M3 1 
halo is not significantly populated with MACHOs. Under this 
assumption, it is expected that the events to be detected to- 
ward the field in and around the M3 1 bulge, where the event 
rate is highest, ar e caused mostly by low-mass stars located 
in the bulge itself (Kerins et al. 2005). We, therefore, choose 
a representative event as the one caused by a lens with a pri- 
mary star mass of m = 0.3 M and distances to the source star 
and lens of D$ = 780 kpc and Di = (780- 1) kpc, respectively. 
Then, the corresponding physical and angular Einstein ring 
radii are r E = 1.56 AU and 9 E = 2.0 fj,as, respectively. The 
assumed timescale is t E = 20 days. 

To see the dependence of the detection efficiency on the 
stellar type of the source star, we test several types of source 
stars with various absolute magnitudes and sizes. The type of 
the source star affects the detection efficiency in two differ- 
ent ways. On one side, high-luminosity of the bright source 
star contributes to the detection efficiency in a positive way 
because the photometric precision improves with the increase 
of the source star brightness. If the source star is too bright, 
on the other hand, it is likely to be a giant star, for which the 
planet ary signal e might be di minished due to the finite-source 
effect ( iBennett & Rhidfl99 61 The stellar types of the tested 
source stars are giant, A5, and F5 main-sequence (MS) stars 
with the absolute magnitudes of Mj = 0.0, 1.73, and 2.86, and 
stellar radii of R+ = 10.0 R@, 1.7 R Q , and 1.3 Rq, respec- 
tively. The corresponding source angular radii normalized by 
the Einstein radius are Pir = 6±/6 E = (R i ,/D s )/e E = 0.03, 0.005, 
and 0.004, respectively. 

Observations and photometry are assumed to be carried out 
as follows. Following the specification of the Liverpool Tele- 
scope, we assume that the instrument can detect 1 photon/s 
for an / = 24.2 star. We also assume that the average see- 
ing is # S ee = l"-0 and the observation is carried out such that 
small-exposure images are combined to make a 30 min ex- 
posure image to obtain a high signal-to-noise ratio while pre- 
venting saturation in the central bulge region. The photom- 
etry is done such that the flux variation is measured at an 
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FIG. 1 . — Contour maps of the detectability of the planetary lensing signal, T>, as a function of the source star position for events caused by planetary systems 
with various star-planet separations, s, and planet/star mass ratios, q, and associated with a giant source star. The source star is assumed to have an absolute 
magnitude M; = 0.0 and a radius = 10.0 Rq . The maps are centered at the position of the primary lens star, and the planet is located on the left. The contours 
are drawn at the levels of 25 = 1.0 (white), 2.0 (yellow), and 3.0 (brown), respectively. The white dotted circle has a radius of uo.th* which represents the threshold 
impact parameter of the source trajectory required for the event detection. All lengths are in units of u® th> whose absolute value is marked in the bottom left panel. 
The closed figure(s) drawn by a blue curve is (are) the caustic(s). The field surface brightness is assumed to be /i = 18.0 mag/arcsec 2 , which is a representative 
value of the M31 bulge region. A sample source trajectory is indicated by the white arrow and the corresponding lightcurve is shown in Fig. 151 

FIG. 2. — Same as Fig.^except that the source star is a A5 main-sequence with an absolute magnitude M/ = 1.73 and a radius /?* = 1.7 Rq. 



aperture that maximizes the signal-to-noise ratio of the mea- 
sured flux variation. In the background-dominated regime 
such as the M31 field, the noise is proportional to the aper- 
ture radius f? ap , i.e. Fb oc n6 2 . On the other hand, assuming 
a gaussian PSF, the measured source flux variation scales as 

F = F S (A- 1) oc f^(9/a 2 PSF )exp{-9 2 /2al SF )de, where a PSF = 
0.425f? see . Therefore, the signal-to-noise ratio scales as S/N = 

F/y7£oc(l/0 a p) Jo i,p (e/aj SP )exp(-0 2 /24 SF ) dO. Then, the 
optimal aperture that maximizes the signal-to-noise ratio is 
#ap = 0.673f? see . With the adoption of this aperture, the frac- 
tion of the source flux within the aperture is F(6 < 9 ap )/F m = 
0.715, where F tot is the flux measured at ap = oo. 

In Figure [00 we present the contour maps of the de- 
tectability of the planetary lensing signal as a function of the 
source star position for events caused by planetary systems 
with various s and q, and involved with source stars of various 
stellar types. The contours are drawn at the levels of T> = 1.0 
(white), 2.0 (yellow), and 3.0 (brown), respectively. We as- 
sume that the planetary signal is firmly detected if T> > 3.0. 
In the map, we present only the region around the 'detection 
zone', which represents the region of the source star position 
where the magnification is higher than a threshold magnifica- 
tion required for the event detection, A±. The threshold mag- 
nification is defined by (A t h — l)F$ = 3\/7b, i-e-, 3<r detection 
of the event. If we define Mo.th as the threshold lens-source 
impact parameter corresponding to Ath, everything of interest 
is contained within the circle with the radius wo.th (marked by 
a white dotted circle in each panel). We, therefore, use Mo.th 
as a scale length instead of the Einstein radius. However, to 
provide the relative size of the detection zone, we mark the 
absolute value of Mo.th in the bottom left panel of each fig- 
ure. The value of the threshold impact parameter decreases as 
either the source star becomes fainter or the background sur- 
face brightness increases. In Figure^ we present the variation 
of the threshold impact parameter as a function of the back- 
ground surface brightness for source stars of various types. 
The maps are constructed for a common surface brightness 
of fi = 18.0 mag/arcmin 2 , which is a representative value of 
the M3 1 bulge region. For the construction of the maps, we 
consider the attenuation of the magnification caused by the 
finite-source effect. 

From the maps, we find that although both the sizes of the 
event detection zone (mo < uqjh) and the region of the plane- 
tary perturbation (D > 3) decrease as the source star becomes 
fainter, the rate of the decrease of the planetary perturbation 
region is smaller than the rate of decrease of the event de- 
tection zone. This is because the planetary perturbation is 
confined to the region around the caustic whose size does not 
depend on the source brightness. As a result, the perturbation 
region occupies a greater fraction of the event detection zone 
with the decrease of the source star brightness. However, this 
does not necessarily imply that the planet detection efficiency 



of events associated with MS stars is higher than that of events 
associated brighter giant stars. This is because despite the 
slow rate of decrease, the perturbation region does decrease, 
and thus picking up the resulting short-duration planetary sig- 
nals for events involved with faint MS stars requires higher 
monitoring frequency. 

Once the maps of the detectability are constructed, we pro- 
duce a large number of light curves of lensing events resulting 
from source trajectories with random orientations and impact 
parameters mo < «o,th ( see example light curves in Figure 0. 
Then, we estimate the detection efficiency as the ratio of the 
number of events with detectable planetary signals to the to- 
tal number of the tested events. We assume that on average 
five combined images with a 30 min exposure are obtained 
daily following the current Angstrom survey. By applying a 
conservative criterion for the detection of the planetary sig- 
nal, we assume that the planet is detected if the signal with 
T> > 3 is detected at least five times during the event. Since 
the monitoring frequency is / = 5 times/day, this implies that 
the planetary signal should last at least 1 day for detection. 

In Figure |6] we present the estimated detection efficiency 
as a function of s and q for events involved with various 
source stars. In Figure0 we also present the variation of the 
efficiency depending on the background surface brightness, 
where the efficiency is estimated by varying /j, but fixing the 
lens parameters as s = 1 .2 and q = 5 x 10~ 3 . From the figures, 
we find the following results. 

1 . For events associated with giant source stars, it will be 
possible to detect planets with masses equivalent to or 
heavier than the Jupiter (q ~ 3 x 10~ 3 ). Although the 
efficiency varies considerably depending on the star- 
planet separation, the average efficiency is ~ 3% for 
events caused by a lens system having a planet with 
a mass ratio q = 5 x 10~ 3 and located in the lensing 
zone. However, it is expected that detecting planets 
with masses less than the Saturn (q ~ 10~ 3 ) would be 
difficult. 

2. The optimal events for the detections of the planetary 
signals are those associated with giant stars. It is ex- 
pected that detecting planets for events associated with 
MS stars would be difficult because of the poor photom- 
etry and the resulting short durations of the planetary 
signals. For example, the duration of the event associ- 
ated with an F-type MS star is fdm- < 2Mo, t hfE ~ 1 day. 
Considering that the planetary perturbation region with 
D > 3 occupies a fraction of the detection zone (as 
shown in Fig. [3}, the planetary signal would be too 
short for detection. Although planets can be detected 
with a non-negligible efficiency (~ 1 %) for events as- 
sociated with A-type MS stars, the number of planet 
detections from these events would be small because of 



4 

Fig. 3.— 



Chung et al. 

Same as Fig.0except that the source star is a F5 main-sequence with an absolute magnitude Mj = 2.86 and a radius = 1 .3 i?Q . 
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FIG. 4. — Threshold lens-source impact parameter, uo.ih. as a function of 
the background surface brightness fj, for source stars of various types. 



the rarity of early-type MS stars projected on the M31 
bulge region. MS-associated events could be detected 
at low surface-brightness regions as shown in FigureQ 
but the event rate toward these fields would be low due 
to the low column density of lens matter along the line 
of sight. 

The efficiency peaks at a certain surface brightness. In 
the region of very high surface brightness, planet detec- 
tion is limited by the poor photometry. In the very low 
surface-brightness regime, on the other hand, the event 
detection zone is significantly larger than the size of the 
planetary deviation region, which is confined around 
caustics. As a result, the efficiency, which is propor- 
tional to the one-dimensional size ratio of the planetary 
deviation region to the detection zone, is low in this re- 
gion. The peak efficiency occurs at [i ~ 20 mag/arcsec 2 
for events involved with giant source stars (see Fig- 
ure 0. 

In most cases, the planetary perturbations of M3 1 pixel- 
lensing events are induced by central caustics. There- 
fore, an observational strategy of focusing on these per- 
turbations would maximize the detections of M3 1 plan- 
ets. An alert system based on real-time survey obser- 
vations combined with prompt follow-up observations 
would do this. 




FIG. 5. — Example pixel-lensing light curves with planetary perturbations. 
The time scale is normalized by the effective time scale t e g that is the time re- 
quired for the source to cross the threshold impact parameter. The individual 
light curves are for events produced by a common planetary lens system with 
s = 1.2 and q = 5 X 10~ 3 , but occurred on different source types of giant, A5, 
and F5 main-sequence stars. The source trajectories of the individual light 
curves are marked in the corresponding panels in Fig. and [3] respec- 
tively. The data points on each light curve are marked under the assumption 
that the monitoring frequency of 5 observations per 24-hour period. The dot- 
ted curve shows the lightcurve expected in the absence of a planet. 



TABLE 1 

Detection Efficiencies Under Various Strategies 



strategy 


source star type 


detection efficiency 


I 


giant 


3% 




A5 MS 


1% 




F5 MS 


0% 


II 


giant 


7% 




A5 MS 


9% 




F5 MS 


1% 


III 


giant 


41% 




A5 MS 


52% 




F5 MS 


66% 



Note. 

n-3 



Average detection efficiencies of detecting a planet with q = 
5 X KT 3 located in the lensing zone of 0.6 < s < 1.6 under three different 
observational strategies. The individual strategies imply a survey mode with 
a monitoring frequency / = 5 times/day (strategy I), a survey mode with 
/ = 10 times/day (strategy II), and a survey mode with / = 5 times/day plus 
follow-up observations by employing a single large telescope (strategy III), 
respectively. 



4. IMPROVING PLANET DETECTION EFFICIENCY 

Considering that planetary perturbations, in many cases, are 
missed from detection due to short durations, a significant im- 
provement in the planet detection efficiency is expected with 
increased monitoring frequency. One way to do this is using 



more telescope time or employing more telescopes ('strategy 
IF). The other way is conducting follow-up observations for 
events detected in the early phase by the survey experiment 
(strategy TIF). In this section, we estimate the efficiencies ex- 
pected under these improved observational strategies. We des- 
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FIG. 6. — The planet detection efficiency of the Angstrom survey as a 
function of the star-planet separation s, and the planet/star mass ratio q, for 
M31 lensing events associated with various source stars. The type of the 
involved source star is marked at the bottom left corner of each panel. The 
field surface brightness is assumed to be = 18.0 mag/arcsec 2 . 



ignate the observational condition of the current pixel-lensing 
survey (with / = 5 times/day) as 'strategy I'. 

We simulate the observations under the strategy II by dou- 
bling the monitoring frequency of the current Angstrom ex- 
periment, i.e. /= 10 times/day. For strategy III, we assume 
survey-mode observations with / = 5 times / day and follow- 
up observations with a single 8m-class telescope. Follow- 
up observations are assumed to be carried out 4 hours after 
the first pixel-lensing event signal is detected by the survey 
observations and the assumed monitoring frequency is / = 
20 times/night. Since a single telescope is employed, follow- 
up observations are assumed to be carried out only during the 
night (8 hrs per day), and thus 24 min per each combined im- 
age. Assuming 20 min per exposure (to allow ~ 4 min for 
readout), the photometric uncertainty of the follow-up obser- 
vation is <78m/c2m ~ (30 min/20 min) I//2 (2m/8m) ~ 31% of 
the survey observation. 

In Figure [8] we present the efficiency expected under the 
two improved observational strategies. In Table [2 we also 
present the average detection efficiencies of detecting a planet 
with q = 5 x 10~ 3 located in the lensing zone under the three 
different observational strategies. From the figure and ta- 
ble, we find that significant improvement in efficiency is ex- 




/J, (mag/ arcsec 2 ) 



FIG. 7. — Variation of the planet detection efficiency depending on the 
background surface brightness /i. The efficiency is estimated for the case of 
a planetary system with the star-planet separation and the planet/star mass 
ratio of s = 1.2 and q — 5 X 10~ 3 , respectively. 



pected, especially from the adoption of the follow-up observa- 
tion strategy. The improvement is more significant for events 
involved with MS stars because the short-duration perturba- 
tions associated with these events are readily detectable with 
the increased monitoring frequency. 

5. CONCLUSION 

We explored the feasibility of detecting planets in M31 
from a high-frequency pixel-lensing survey using a global 
network of 2m-class telescopes. For this evaluation, we esti- 
mated the efficiency of detecting planetary signals for events 
induced by planetary systems with various planet/star mass 
ratios and star-planet separations, associated with source stars 
of various types, and occurring toward fields with various sur- 
face brightness. From the dependence of the detection effi- 
ciency on these parameters, we found that 3cr detection of the 
signals produced by giant planets located in the lensing zone 
with masses equivalent to or heavier than the Jupiter would 
be possible with detection efficiencies ~ 3% if the event is 
associated with giant stars. A dramatic improvement of the 
efficiency is expected if follow-up observations based on real- 
time survey observations become possible. 
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FIG. 8. — Same as Fig. [6] but the efficiency is estimated under improved observational setups. The left and right panels are the efficiencies estimated under the 
observational strategy II and III, respectively. See the details of the observational strategies in § 4. 



and do not necessarily reflect the view of the NSF. 
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